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ABSTRACT 

Aims. We study chromospheric emission to understand the temperature stratification in the solar 
chromosphere. 

Methods. We observed the intensity profile of the CallH line in a quiet Sun region close to the 
disk center at the German Vacuum Tower Telescope. We analyze over 10^ line profiles from inter- 
network regions. For comparison with the observed profiles, we synthesize spectra for a variety of 
model atmospheres with a non local thermodynamic equilibrium (NLTE) radiative transfer code. 
Results. A fraction of about 25 % of the observed Ca ll H line profiles do not show a measurable emission 
peak in H2v and H2t wavelength bands (reversal- free) . All of the chosen model atmospheres with a 
temperature rise fail to reproduce such profiles. On the other hand, the synthetic calcium profile of a 
model atmosphere that has a monotonic decline of the temperature with height shows a reversal-free 
profile that has much lower intensities than any observed line profile. 

Conclusions. The observed reversal-free profiles indicate the existence of cool patches in the interior 
of chromospheric network cells, at least for short time intervals. Our finding is not only in confiict 
with a full-time hot chromosphere, but also with a very cool chromosphere as found in some dynamic 
simulations. 
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1. Introduction 

Solar model atmospheres are built to reproduce the 
observed intensity spectrum of the Sun in contin- 
uum windows and se lected lines (e.g., th e series 
of models initiated bv lVernazza et anil98ll ). These 
models rely on temporally and spatially averaged 
spectra. Hence, they provide information about the 
average properties of the physical parame ters in the 
solar atmosphere ()Linskv fc Avrettlll970[ ). At high 
spatial, spectral, and temporal resolution, differ- 
ences from the average profiles arise, leading to dif- 
ferent interpretations of the atmospheric stratifica- 
tion. In this context, there is no general agreement 
on the temperature stratification in the solar chro- 
mosphere. This is closely related to the fundamen- 
tal proble m of the heating mechan i sms in the chro- 
mosphere. ICarlsson fc SteinI (|1994 1X9971) proposed 
that apart from short-lived heating episodes in the 
chromosphere its time-averaged temperature strat- 
ification is monotonically decr easing w ith geometri- 
cal height. As noted bv .Carlsson et al.l (|l997), there 
are disagreements between this model and SUMER 
data. Moreover , the cool model was criticized by 
iKalkofen et all (|1999f ) who argued that the chro- 



mosphere is always hot and can never reach such a 
low activity state a s proposed by Carlsson fc Stein. 
On the other hand, lAvred ((2002) challenged the ar- 
guments of Kalkofen et al. based on observations of 
CO lines at the solar limb. The results of the CO 
observations were supported by three-dimensional 
MHD simulations that suggested very cool struc- 
tures in the upper atmosphere and implied a ther- 
mally bifurcated medium in non-ma gnetic regions 
(jWedemever Bohmet aLll2004 l2005l) . 

Recently, lAvrettl (|2007 *) suggested that it is 
possible to reproduce a wide range of continuum 
and line intensities with a temperature stratifi- 
cati on within 400 K of the semi-empirical mod- 
els ()Fontenla et al.l l2006l and references therein) . 
He concluded that "...these results appear to con- 
flict with dynamical models that predict time vari- 
ation of 1000 K or more in the chromosphere...". 

In this paper, we comment on both full-time hot 
and dynamic models, based on the lowest intensity 
profiles we observed. We present observations and 
data reduction in Sect. 2. Then, we try to reproduce 
the observed reversal- free calcium profiles (i.e., pro- 
files without emission peaks) by synthesizing line 
profiles from a variety of model atmospheres using 
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Fig. 1. Maps of the observed area. Top: H-index in pm; middle: profile type. It is coded as follows: 0: 
reversal- free; 1: norma l; 2: asymmetric; 3: other types. Bottom: Vtot, which is a measure of the net vertical 
magnetic flux density (jLites et al.lll999l ) in arbitrary units. The a;-axis shows both scanning direction and 
time (width of a single map is eight arcsec). The y-axis is the slit direction. Both axes are in arcsec. Only 
17 out of 40 scans are shown here. The time duration per scan was about 100 s. 
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a non local thermodynamic equilibrium (NLTE) ra- 
diative transfer code (Sects. 3 and 4). The discus- 
sion and conclusion appear in Sects. 5 and 6, re- 
spectively. 

2. Observations 

A quiet Sun area close to disk center (cos 6* = 0.99) 
was observed at the German Vacuum Tower 
Telescope (VTT) in Tenerife, July 07, 2006. The 
good and stable seeing condition during the ob- 
servation and the Kiepen heucr Adaptive Optics 
System (jvon der Liihe et a l. 2003) provided high 
spatial resolution in the calcium profiles, leading to 
a collection of different structures observed in the 
quiet Sun. We used an integration time of 4.8 s and 
achieved a spatial resolution o f 1 arcsec. Details of 
the observation appeared in Rezaei et al.l (|2007bl ). 

The CallH intensity profile and Stokes vec- 
tor profiles of the visible neutral iron lines at 
630.15 nm, and 630.25 nm were observed with 
the blue and red channels of the PO larimctric 
LIttrow Sp e ctrog raph (POLIS: lSchmidt et al. 2003t 
iBeck et al.l 120051 ). POLIS was designed to pro- 
vide co-temporal and co-spatial measurements 
of the magnetic field in the photosphere and 
the CaiiH intensity profile. We normalized the 
average calcium profile in the line wing at 
396.49 nm to the Fourier Transform Spectrograph 
profile (jStenflo et al.lll98l ) and applied the calibra- 
tion factor to all profiles. Details of the calibration 
procedure fo r the ca lcium profiles are explained in 
iRezaei et al.l (|2007al) . The overall properties of the 
observed time-series are shown in Fig. [1] concate- 
nated to a single map. The H-indcx, which is the 
integral of the normalized calcium profile within 

0. 1 nm around the core, i s shown in the to p panel 
fsee lCram fc Damelll983l: iLites et al.lll993l for def- 
initions of the H-index, V and R band intensi- 
ties). The middle panel shows the profile type. It is 
coded as follows: [0] reversal-free, [1] normal (with 
two emission peaks), [2] asymmetric (with only one 
emission peak), and [3] other types. The bottom 
panel of Fig. [T] shows V tpt, which is a m easure of 
the signed magnetic flux (jLites et al.lll999|) . It is im- 
mediately seen that reversal-free profiles are absent 
in the vicinity of network elements {y ~ 20 arcsec 
in Fig.IH). 

The top row of Fig. [2] shows three profiles be- 
longing to different activity levels in the quiet Sun. 
The left profile has two red emission peaks, the 
middle profile shows virtually no violet emission 
peak and the right profile shows a strong violet 
and a weak red emission peak {V/R > 2.1). The 
gray profile is the average of a few thousand pro- 
files (cell interior and boundary) with an H-index of 
10.1pm. An alternative way to measure the asym- 
metry of the two emission peaks is dividing the 
pea k values rather tha n band integrated intensi- 
ties (jRezaei et al.ll2007al) . In very asymmetric cases, 

1. e., where there is only a large violet peak, we can 
easily achieve a large V/R ratio with band defini- 



tions whereas the peak ratio is undefined. Double 
reversal calcium profiles in old observations with 
lower spatial resolution do not show extreme val- 
ues for the ratio of the peaks. It requires both 
high spatial and temporal resolution and large sig- 
nal to noise ratio. The maximum and minimum 
V/R (peak) ratio in our data set are 2.2 and 0.5 
(top right and middle panels. Fig. [5]). The impor- 
tance of using peak values instead of bands is that 
the formation height of different wavelengths in the 
band can differ by more than one Mm (see Fig. [5] 
and the left panel of Fig. [6]). The bottom panels 
of Fig. [5] show calcium profiles observed in abso- 
lutely quiet Sun regions. There is no emission peak 
or bulge in the core of these profiles. In the bot- 
tom row of Fig. [21 the left profile has the same 
wing intensity as the average profile (gray), while 
the middle and right profiles show lower and higher 
wing intensities than the average profile, respec- 
tively. Note that the reversal-free profiles usually 
have a V/R ratio very close to one, i.e., they are 
symmetric (check values in the bottom panels of 
Fig. [2). These reve rsal-free profiles were ca l led ab - 
sorption profiles bv I Grossman- Doerth et al](|1974D . 
They rightly concluded that the source function of 
the calcium line then shows negligible enhancement 
in the atmospheric layers where the line forms. An 
emission bulge exists in the minimum profile of 
ICram fc Daini ([1981) (their Fig. O while this is 
not the case in the bottom panel, Fig. O This dif- 
ference means that our reversal-free profiles present 
a new minimum state for the observed Ca ii H pro- 
files. Our observations have better spatial an d tem- 
poral resolutions than lCram fc Dami ()1983f ). so we 
expect to find more extreme cases than these au- 
thors. 



Statistics: There are many references on statis- 
tical properties of the CaiiH line profiles (e.g., 
ICram fc Damel[l983D , but our statistics on reversal- 
free profiles are the first. We exclude the range 
2/ = 15 — 25 arcsec in our maps because that clearly 
contains magnetic network. Then we perform a sta- 
tistical analysis to check how often the reversal-free 
profiles occur. The inter-network sample has some 
129 000 profiles. We check the existence of the emis- 
sion peaks by analyzing the profile shape. Among 
the profiles that have no emission peaks, some show 
a strong H-index as in the top panel of Fig. [U 
but others, like the bottom panels of Fig. [21 show 
a rather normal wing intensity and a dark core. 
We find that 25 % of all profiles do not show any 
emission peak: 16% have an H-index larger than 
9 pm, 8 % between 8 and 9 pm, and 1 % below 8 pm. 
Examples of calcium profiles with low H-index and 
no emission peak are shown in the lower panels of 
Fig. [2l A majority of these 25 % show a tiny bulge 
at the position of the emission peaks. 

It is worth noting that the H-index alone is not a 
reliable parameter for identifying the reversal-free 
profiles, as is demonstrated in Fig. [31 The upper 




Fig. 2. Examples of single calcium profiles with drastically different properties observed in the quiet 
Sun. The gray profile is the average of a few thousands profiles (cell interior and boundary). Top left: a 
profile with three emission peaks. Top middle: a very asymmetric profile with only a red emission peak. 
Top right: a very asymmetric profile with two emission peaks. Bottom: three reversal-free profiles. Note 
the level of the wing intensity with respect to the wing of the average profile. The dashed lines show a 
distance of ± 0.05 nm from the average calcium core, indicated by a vertical dotted line. The bottom left 
profile is an average of 59 reversal-free profiles while the other panels show single profiles. Except for the 
top middle and right panels, the V/R ratio was calculated using bands. 



profile shows no emission peaks, but the wing in- 
tensities suggest a hotter upper photosphere and a 
lower photospheric temperature gradient; the lower 
profile has two clear emission peaks, but the wing 
intensities demand a cooler upper photosphere. 
Consequently, the H-index of the bottom profile is 
less than half of the upper profile. This is in accor- 
dance with the fact that all three panels of Fig. [2] 
(bottom row) show reversal-free profiles with an H- 
index larger than 8 pm. 

Besides the 25 % fraction of reversal-free profiles 
in our sample, 46% of all profiles show a double 
reversal, with strong tendency to have a brighter 
violet peak. The remaining 29 % of all inter-network 
profiles show just one emission peak, again with a 
majority of strong vio let emission p eaks, similar to 
previous studies (e.g.. lCram fc Dame 1983 ). 



3. Synthetic calcium profiles 

To reproduce the reversal-free profiles we compute 
synthetic Ca II H profiles by means of t he NLTE ra- 
diativ e transfer code RH, developed bv lUitenbroekl 
(|200lh . The atomic model is the standard 5-level 
plus continuum model that is included with the 
RH code; apart from minor atomic data updates, 
thi s model is essentially equal to the one used 
by lUitenbroekl fl989„), w hich in turn dates back 
to IShine fc Linskvl (119741) and references therein. 
As shown bv lUitenbroek (|1989D . partial frequency 
redistribution of the photons in both CaiiH&K 



needs to be taken into account, whereas assuming 
complete redistribution in the lines of the infrared 
triplet has no noticeable influence on the profiles of 
H&K. Partial redistribution makes the line source 
functions wavelength dependent and significantly 
reduces the strength of the H2 and K2 emission 
w.r.t. the common complete redistribution case. 
Note that since we only deal with static atmosphere 
models, it suffices to use angle-averag ed partial fre- 
quenc y redistribution. Furthermore, ISolanki et al.l 
(| 199 If) found that the width of the K2 (and im- 
phcitly H2) peaks correlates with the microturbu- 
lent velocities in the chromosphere. Zero micro- 
turbulence produces stronger and narrower peaks 
whereas the microturbulence values given with 
semi-empirical ID solar models typically result in 
shallower, wider peaks. Microturbulence is a pa- 
rameter that is used in particular in ID radiative 
transfer modeling to account for the influence of 
spatially unresolved velocities on the line profiles. 
For 2D and 3D dynamic models it is generally not 
needed. 

Since full 2D or 3D (dynamic) radiative trans- 
fer modeling is beyond current computational re- 
sources, we use static ID plane-parallel models to 
synthesize CaiiH line profiles, even though the 
structures that we observe are so small that we ex- 
pect lateral radiative transfer to play a role. In ad- 
dition, we use the microturbulence values as speci- 
fied with the semi-empirical models we employ, even 
though we do not know anything about the un- 
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Fig. 3. These two examples show that the inte- 
grated intensity close to the core, either H-index 
or a filter, is not a reliable criterion to separate 
reversal-free profiles from profiles with emission 
peaks. Positions (in Fig. [T]) of the upper and lower 
profile are (18.7, 30.5) and (32.6, 30.2), respectively. 



resolved small-scale velocities in the observed fea- 
tures. Synthetic profiles from semi-empirical ID 
quiet-Sun models with non-zero microturbulence 
tend to have H2-peaks that are generally stronger 
than observed; using zero microturbulence, as one 
would do if the structure to be modeled were com- 
pletely resolved, would only widen this discrepancy. 
Since we cannot be sure that the structure is com- 
pletely resolved in the observations, and also be- 
cause we expect that radiative interaction with the 
environment will lead to broadening of the H2- 
peaks, we decided not to change the microturbu- 
lent velocities. Note that the latter is not an obser- 
vational effect, but solely due to scattering of line 
photons in the solar atmosphere. 

As model atmospheres we use four semi- 
empirical models: FALC (quie t Sun) and 
FALA (inte r-network) from iFontenla et al.l 
(|l999l l2006t): the Holw e ger-M iiller atmo- 
sphere ( Holweger fc Miilleil [TItI), which is 
very similar to a theoretical radiative equilib- 
rium model, but not constructed to strictly 
satisfy this condition; and the new C6 quiet 



Sun model a.tmosphere, prov i ded in advance of 
publication (jAvrett fc Loeseil l2008l where this 
model is now called C7). The left panel of Fig. |4] 
shows the temperature stratification of these four 
model atmospheres. All of the models, except the 
Holweger-Miiller model, show a chromospheric 
temperature rise at different heights in the at- 
mosphere. The Holweger-Miiller model has a 
monotonically-decreasing temperature stratifica- 
tion with an artificial extension with constant slope 
(almost isothermal). The extension is required to 
include the formation height of the Ca 11 resonance 
lines. The synthesized calcium profiles are shown 
in the right panel of Fig. [H For comparison, the 
reversal-free profile from the bottom left panel of 
Fig. [21 is overplotted in gray. 

For all synthesized profiles, except the one de- 
rived from the Holweger-Miiller model, there are 
well-defined emission peaks at the and i?2r 
wavelengths. All the synthetic profiles are signifi- 
cantly different from the observed reversal-free pro- 
files. It might suggest that, e.g., the temperature 
rise in the chromosphere is either shifted to higher 
layers, or the temperature gradient is not as steep 
as proposed by, for example, the C6 model. For a 
very inhomogeneous model atmosphere, the region 
with high emission may be very "patchy" . 

The new chromospheric model C6, like FALC, 
includes a chromospheric temperature rise. If one 
considers a velocity stratification in the atmo- 
sphere, it will be possible to reproduce a vari- 
ety o f asymmetric calcium profiles (e.g., ' HeaslevI 
II975D . However, it is not trivial to diminish both 
emission peaks at the same time to reproduce a 
reversal-free profile. As seen in Fig. [4l none of 
the applied models with a chromospheric temper- 
ature rise produces a reversal-free calcium profile. 
While the C6, FALA, and FALC models look too 
hot, with strong emission peaks, the Holweger- 
Miiller model fails due to its v ery de ep calci um 
core. The d yna mical models of Carlsson fc S teinI 
(|1994I1997D and lRammacher fc Cuntzl (|200a) have 
episodes (significantly) cooler than the Holweger- 
Miiller mo del with correspondingly lower line c ore 
intensities (|Uitenbroekll2002t lRammacherj |2005') . It 
should be remembered, however, that the chromo- 
spheric radiative transfer in those simulations is not 
"consistent" : important cooling agents are not rep- 
resented at all and deviations from various equilib- 
rium states are ignored. 

4. Modified model atmospheres 

Given the large deviations between all synthetic 
profiles of the Ca 11 H line and the observed reversal- 
free profiles, it is tempting to try to obtain a closer 
match to the observed profiles by making changes 
to one of the model atmospheres, specifically to 
the FALC model. This procedure is by no means 
unique and given that we employ static models it 
is not clear whether we can reproduce the observed 
profiles that result from the dynamic and highly- 
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Fig. 4. Left: Stratification of the temperature in the FALA, FALC, Holweger-Miiller, and C6 model 
atmospheres. Right: Resulting calcium profiles. Except for the Holweger-Miiller model atmosphere, which 
does not show a temperature rise, all model atmospheres result in two emission peaks. We set the nearby 
pseudo-continuum of the FALC profile to unity and use this normalization factor for all synthetic profiles. 
The observed reversal- free profile in the bottom left panel of Fig. [5] is shown in gray. All synthetic profiles 
convolved to the spectral resolution of the observed spectra; scattered light removed from the observed 
profiles. 
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Fig. 5. Intensity contribution function for the in- 
ner part of the Can H line computed for the FALC 
model. The black line is the calcium profile of the 
model. 

structured solar atmosphere. We could end up ad- 
justing all atmospheric parameters and still not get 
a satisfactory match. Instead, we choose to adjust 
only the FALC temperature in a very schematic 
way to get at least an indication which changes 
might be needed. We stress, however, that these 
modified models are no longer consistent since only 
the temperature is changed and all other param- 
eters are left intact. They may not reproduce ob- 
served profiles of other lines. From the NLTE ra- 
diative transfer computations we know that the 
monochromatic source function follows the local 
temperature fairly well throughout most of the 
CaiiH line, even though it is a strong, scattering 
line. This means that a naive change of the temper- 
ature stratification at some height translates into 
a change in emergent intensity for those parts of 



the line that form at that particular height. The 
intensity contribution function for the Ca ii H line 
(Fig. [5|) indicates that the formation range for any 
given wavelength is very limited, which means that 
temperature corrections at a given height will lead 
to predictable line profile changes. This procedure 
should work particularly well for the line wings up 
to the Hi minima, where the source function is 
very close to LTE, and for the very core, whose 
source function drops way below the Planck func- 
tion in a predictable way. However, even though 
the formation range at any wavelength within the 
H2 peaks is narrow, the wide range of formation 
heights spanned by these peaks as a whole necessi- 
tates temperature changes over a wide height range 
to suppress the emission peaks completely. In ad- 
dition, we need a rather high temperature in the 
upper chromosphere to obtain a reasonable core in- 
tensity; for that reason we use the FALC model as 
starting point. This means that we may only change 
the temperature structure in the lower and middle 
chromosphere: we created models with the temper- 
ature minimum region isothermally extended up- 
ward to column mass values of 10~'^ " (ISO_30), 
10-3-5 (ISO_35), and IQ-^o (ISO_40), followed by 
a linear r(logm) relation up to a column mass of 
m — 2 - 10~^ g/cm^, above which the original FALC 
temperatures are retained (Fig. [6l left panel). The 
lower temperature in the height range about 600- 
1900 km weakens the emission peaks, while the tem- 
perature rise in the higher layers prevents a very 
deep calcium core. The profiles (Fig. [51 right panel) 
suggest that a temperature rise shifted to larger 
heights is a step toward producing a reversal-free 
profile. 
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Fig. 6. Left: Stratification of the temperature in the FALC model atmosphere and three modified versions 
of that. The thin curves represent the monochromatic source functions of these models at AA = — 0.014 nm 
from line center. The diamonds (from left to right on each curve) represent the location of = 1 
for AA — —0.016, AA — —0.014, and AA — — 0.012 nm. The source function is similar for all three 
wavelengths, so that we restrict to just one source function curve per model. Right: Resulting calcium 
profiles of the mentioned models convolved to the spectral resolution of POLIS. The ISO_30 profile is 
much closer to the observed reversal-free profile than the FALC model. 



5. Discussion 

We find that a quarter of all observed profiles are 
reversal-free profiles. We observed reversal-free pro- 
files of the CallH line with an exposure time of 
about 4.8 s. However, they may get lost at longer 
exposure times. The maximum spatial extension 
of these reversal-free profiles in our map is about 
5 arcsec; they mostly appear in smaller structures 
(Fig. [U middle panel). The coherency of the posi- 
tion of the reversal-free profiles shows that they are 
not an artifact due to effects caused in the Earth 
atmosphere. The observed reversal-free CaiiH pro- 
files (Fig. [21 lower panels) suggest that there may 
be some cool patches in the chr omosphere for 
short periods of time. In contrast, ICarlss on et al.l 
(|1997[ ) tried to find reversal-free profiles in chromo- 
spheric spectral lines observed with SUMER (46.5- 
161 nm). On the basis of a few lines they concluded 
that "all chromospheric lines show emission above 
the continuum everywhere, all the time" . These 
lines cover a large height range above the classical 
temperature minimum and should also reveal the 
cool patches seen in Ca ii H line profiles if observed 
at the proper spatial and temporal resolution. 

We propose that at a spatial resolution of one 
arcsecond, the average temperature stratification in 
the solar atmosphere is probably hotter than the 
Holweger-Miiller model atmosphere, but still does 
not show a permanent temperature rise at all spa- 
tial positions. Moreover, the residual core intensity 
of these profiles is about 0.051c. The Holweger- 
Miiller model, as well as compa r able cool mod- 
els (like m odel Q of Solanki et al.l (jl99l l) or model 
COOLC of lAvres et al. (1986)) and the cool phases 
of the dyna mic models result in a very dark cal- 
cium core (|Uitenbroekl I2002D . far below the ob- 
served reversal-free profiles. In these models, we 



have a longer cool p hase and a shorter hot phase 
in th e chromosphere ijRuttenl [19981 1 Judge fc Peter! 
|1998[ ). Statistics of our inter-network sample indi- 
cates that some 75 % of all profiles show one or 
more clear emission peaks. This appears to be in 
conflict with these models. But we recall that even 
if at any given point in the atmosphere the tem- 
perature is low most of the time, the CallH line 
may still show emission peaks most of the time be- 
cause it is an integral property with contributions 
from a large height range. We speculate that this 
is even more important in 3D models, where lat- 
eral interaction with the ubiquitous shocks plays a 
role as well. Unfortunately, so far no one has at- 
tempted to compute Ca ii H profiles simultaneously 
with the 3D (M)HD simulations. Even worse, com- 
puting Ca II H a posteriori for a single 3D snapshot 
already presents an unsolved task. 

On the other hand, our statistical studies of the 
inter-network calcium profiles indicate that about 
25 % of the profiles do not show any emission peak 
(but might show a bulge). These reversal- free pro- 
files (bottom panels, Fig. ^ are not reproduced 
by any of the hydrostatic model atmospheres with 
a chromospheric temperat ure rise, like in the new 
C6 model. This challenges [Avrett fc Loeseil (|2008[ ). 
who claim that they can reproduce all contin- 
uum and line intensities with models that have 
a temperature variation of at most 400 K. In ad- 
dition, the observed reversal-free profiles cannot 
be derived f rom the dynamica l models during the 
cool phase (lUitenbroekl l2002l Fig. 6). Similarly, 
lAvrettl (|1985 ) tried to modify the temperature 
stratification of VALA, VALC and VALF models 
to reproduce minimu m and maximum profiles of 
ICram fc Damel (|1983() . While we kept the temper- 
ature minimum constant and extended it isother- 
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mally, they introduced a hotter temperature mini- 
mum and decreased the temperature gradient (their 
Fig. 17). Our model ISO_30 and the correspond- 
ing profile (Fig. [5]) better fits the observed reversal- 
free profile, which has lower intensity at the emis- 
sion peak wave l ength s than the minimum profile of 
iCram & Dai^^ (1983"). 

[Fontcnla ct al. (2007) presented a compromise 
solution for the cool patches. They proposed that 
cool patches may occur in a limited height range, 
such that they contribute substantially to the tem- 
perature sensitive molecules, but affect other spec- 
tral lines very little. This a rgument may apply to 
high resolution filtergrams (jWoger et al.|[200l . If 
we had cool temperat ures in a very limited height 
range as proposed by iFontenla et al.1 (|2007( ). and 
a hot atmosphere elsewhere, we would see emis- 
sion features somewhere in each Ca ii H spectrum, 
which is not the case. Therefore, although we have 
strong indications of cool patches in the solar chro- 
mosphere, they a re more extend e d (in height) than 
was proposed by Fontcnl a et all (j2007[ ). 

p3alasubramaniam (2 00l| ) presented multiple re- 
versal calcium profiles and interpreted them as a 
signature of magnetic flux emergence. We observe 
some calcium profiles with more than two emission 
peaks (Fig. [21 upper left panel) in the quiet Sun, 
without signature of magnetic flux e mergence in 
the p olarimetric channel of POLIS (see iBeck et aTl 
l2007l for a description of this channel). Therefore, 
his explanation is not unique. 

6. Conclusion 

We find that a quarter of the observed calcium 
profiles in the quiet Sun inter-network show a 
reversal-free pattern. We interpret these profiles 
as strong indication for a temperature stratifica- 
tion cooler than that of the average quiet Sun. 
Although they do not imply a very cool atmo- 
sphere, like the one based on CO line observa- 
tions dAv res 20 0j) an d numerical simulation s (e.g ., 
ICarlsson et all 119971 : iRammacher fc Cuntd 120051 ). 
it should be clea r ly coo ler than was suggested by 
lAvrett &: Loeseil (|2008[ ). The reversal- free profiles 
we observe, at a spatial resolution of 1 arcsec and 
a temporal resolution of 5 s, show a residual inten- 
sity at the core of about 0.051c. Hence, the profiles 
suggest that the cool components are cooler than 
the FALA, FALC, and C6 models in the sense that 
their corresponding profiles do not have any emis- 
sion peak, but hotter than the extended Holweger- 
Miiller model or the cool phase in the dynamic mod- 
els. 
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